---------------------------------------------------------------------------------------------------------------------
       log:  D:/data/23211 Todesursachenstatistik/Nutzungen/1813-2011 ifo Muenchen - Knoche//Analysen/20121004/201210
> 04_tus-analyse.log
  log type:  text
 opened on:   4 Oct 2012, 14:26:16

. 
. *** Titel des Projekts: Analyse des Einflusses von Religionszugehoerigkeit auf natuerliche und gewaltsame Todesursa
> chen 
. *** Datengrundlage:     <Todesursachenstatistik> 
. *** 
. *** Dateiname des Programmcodes: tus-analyse.do
. *** erstellt: 06.08.2012
. *** von: Sascha O. Becker und Ludger Woessmann
. *** E-Mail: sascha.becker@gmx.de
. *** Tel. 0044 7878 7578 11
. *** 
. *** Dateiname des Output-Files: tus-analyse.log 
. *** 
. *** 
. *** Grundriss des Programms: Untersuchung von Religionszugehoerigkeit und Selbstmord
. *** 
. *** Verwendete Variablen:
. *** Originalvariablen:  
. *** ef304             Wohnort der/des Verstorbenen
. *** ef305u3           Sterbedatum - Jahr
. *** ef306             Geschlecht
. *** ef310             Familienstand
. *** ef313             Religionszugehoerigkeit
. *** ef314             Staatsangehoerigkeit
. *** ef316             Todesursache (ICD-Schluessel)
. *** ef317             Todesursache (Ursache der Schaedigung)
. *** ef320u3           Sterbealter - in Jahren
. ***
. *** Neu angelegte Variablen:    Stetige Variablen:
. *** ags               Amtlicher Gemeindeschluessel Wohnort der/des Verstorbenen numerisch statt string
. *** age               Sterbealter
. *** age2              Sterbealter^2 / 100
. *** Dummy-Variablen:
. *** male              maennlich
. *** famstand_unknown  Familienstand unbekannt 
. *** famstand_single   Familienstand Single
. *** famstand_married  Familienstand verheiratet
. *** famstand_widowed  Familienstand verwitwet 
. *** famstand_divorced Familienstand geschieden
. *** rel_unknown       Religionszugehoerigkeit unbekannt
. *** rel_protestant    Religionszugehoerigkeit protestantisch
. *** rel_catholic      Religionszugehoerigkeit katholisch oder orthodox
. *** rel_otherchrist   Religionszugehoerigkeit sonstige christliche
. *** rel_jewish        Religionszugehoerigkeit juedisch
. *** rel_islam         Religionszugehoerigkeit muslimisch
. *** rel_other         Religionszugehoerigkeit sonstige
. *** rel_none          Religionszugehoerigkeit gemeinschaftslos
. *** rel_pro_no        Religionszugehoerigkeit protestantisch oder gemeinschaftslos
. *** german            Staatsangehoerigkeit Deutsch
. *** suicide           Selbstmord-Dummy
. ***
. *** auf Kreisebene (ags) aggregierte Variablen, die als Kontrollvariablen herangezogen werden
. *** n_protestant_ags: Anzahl Protestanten an den Toten im Landkreis ags
. *** n_catholic_ags: Anzahl Katholiken an den Toten im Landkreis ags
. *** prot_in_protcath_ags: Anteil Protestanten an den verstorbenen Katholiken und Protestanten des Landkreises ags
. ***
. *** Variablen aus der VZ1987 auf Landkreisebene:
. *** log_pop_tot_res_t             log(Bevoelkerungszahl im Landkreis)
. *** prot_in_protcath              Protestantenanteil (an Christen)
. *** share_pop_under3_res_t        Bevoelkerungsanteil Alter 0-2
. *** share_pop_3_5_res_t           Bevoelkerungsanteil Alter 3-4
. ***    etc etc
. *** share_pop_70_75_res_t         Bevoelkerungsanteil Alter 70-74
. *** share_pop_75andmore_res_t     Bevoelkerungsanteil Alter 75 plus
. *** share_pop_tot_res_single_t    Anteil Singles an Gesamtbevoelkerung
. *** share_pop_tot_res_married_t   Anteil Verheiratete an Gesamtbevoelkerung
. *** share_pop_tot_res_wid_t       Anteil Witwen/Witwer an Gesamtbevoelkerung
. *** share_pop_tot_res_div_t       Anteil Geschiedene an Gesamtbevoelkerung
. *** share_foreigners              Auslaenderanteil an Gesamtbevoelkerung
. *** share_pop_res_wrkmninc_1st_t  Anteil Gesamtbevoelkerung, die arbeitet in Land- und Forstwirtschaft, Fischerei
. *** share_pop_res_wrkmninc_2nd_t  Anteil Gesamtbevoelkerung, die arbeitet in Produzierendes Gewerbe
. *** share_pop_res_wrkmninc_3rd_t  Anteil Gesamtbevoelkerung, die arbeitet in Handel, Verkehr und Nachrichtenuebermi
> ttlung
. *** share_pop_res_wrkmninc_oth_t  Anteil Gesamtbevoelkerung, die arbeitet in uebrige Wirtschaftsbereiche
. *** share_pop_res_incsub_t        Anteil Gesamtbevoelkerung, die lebt von Arbeitslosengeld/-hilfe, Rente, Pension
. *** share_pop_res_incpar_t        Anteil Gesamtbevoelkerung, die lebt von Unterhalt durch Eltern, Ehegatten usw.
. *** share_school_prim             Anteil Gesamtbevoelkerung mit Volksschule, Hauptschule
. *** share_school_lowsec           Anteil Gesamtbevoelkerung mit Realschule oder gleichwertiger Abschluss 
. *** share_school_uppsec           Anteil Gesamtbevoelkerung mit Hochschul-/ Fachhochschulreife
. *** share_voc_fach                Anteil Gesamtbevoelkerung mit Berufsfach-/ Fachschule
. *** share_voc_uni                 Anteil Gesamtbevoelkerung mit Hochschule / Fachhochschule
. *** dist_wittenberg               Entfernung nach Wittenberg
. 
. **** Lege Arbeitsverzeichnis an
. *capture mkdir "$root/Analysen/20120619"
. 
. *** Erlaube Durchlauf des Programms
. set more off

. 
. /*******************/
. /* 1992-1997: ICD9 */
. /*******************/
. 
. *** Datensatz einlesen
. use "$root/KDFV/tus080_1992.dta", clear

. 
. *** Generiere Amtlichen Gemeindeschluessel Wohnort der/des Verstorbenen numerisch statt string
. gen ags=real(ef304)

. 
. *** Generiere Dummyvariable fuer Maenner
. gen male=(ef306==1)

. 
. *** Generiere Dummyvariablen fuer Familienstand
. gen famstand_unknown  = (ef310==0 | ef310==9) 

. gen famstand_single   = (ef310==1)

. gen famstand_married  = (ef310==2) /* reference category */

. gen famstand_widowed  = (ef310==3)

. gen famstand_divorced = (ef310==4)

. 
. *** Generiere Dummyvariablen fuer Religionszugehoerigkeit
. gen rel_unknown = (ef313==0) 

. gen rel_protestant = (ef313==1 | ef313==2)

. gen rel_catholic = (ef313==3 | ef313==4 | ef313==5) /* reference category */

. gen rel_otherchrist = (ef313==6)

. gen rel_jewish = (ef313==7)

. gen rel_islam = (ef313==8)

. gen rel_other = (ef313==9 | ef313==10)

. gen rel_none  = (ef313==11)

. gen rel_pro_no = rel_protestant + rel_none

. 
. *** Protestantenanteil an "Protestanten plus Katholiken" auf Landkreisebene aus TUS
. bysort ags: egen n_protestant_ags = sum(rel_protestant)

. bysort ags: egen n_catholic_ags = sum(rel_catholic)

. bysort ags: gen prot_in_protcath_ags = n_protestant_ags/(n_catholic_ags+n_protestant_ags)

. 
. *** Deutsche vs Auslaender
. gen german = (ef314==0)

. 
. *** Generiere Selbstmord-Dummy
. *** Selbstmord-Codes in ICD9 sind Codes 950 bis 959 in den FDZ-Daten, aber E950 bis E959 in ICD9-Doku, d.h. 10 Code
> s
. gen suicide=0

. foreach num of numlist 950/959 {
  2.   replace suicide=1 if ef317=="`num'"
  3. }
(* real changes made)
(* real changes made)
(* real changes made)
(* real changes made)
(* real changes made)
(* real changes made)
(* real changes made)
(* real changes made)
(* real changes made)
(* real changes made)

. 
. *** Generiere alter^2 auf Basis von Alter: Sterbealter - in Jahren
. rename ef320u3 age

. gen age2 = age * age / 100

. 
. *** Sortiere nach AGS
. sort ags

. 
. *** Verknuepfe mit Kreisdaten der VZ 1987
. merge ags using "$root/KDFV/VZ1987.dta"
variable ags does not uniquely identify observations in the master data

. 
. *** Benenne Variable Protestantenanteil um
. ren share_prot_in_protcath prot_in_protcath

. 
. *** Betrachte Korrelation zwischen Protestantenanteil in TUS und Protestantenanteil in VZ1987
. pwcorr prot_in_protcath prot_in_protcath_ags, sig

             | prot_i~h prot_i~s
-------------+------------------
prot_in_pr~h |   1.0000 
             |
             |
prot_in_pr~s |   0.9958   1.0000 
             |   0.0000
             |

. 
. *** Generiere log(Bevoelkerungszahl) im Landkreis
. gen log_pop_tot_res_t = log(pop_tot_res_t)
(* missing values generated)

. 
. *** Tabelliere Religionszugehoerigkeit der Verstorben insgesamt, sowie fuer die Selbstmoerder
. *tab ef313
. *tab ef313 if suicide==1
. 
. 
. 
. *** Regressionen: was "erklaert" Selbstmord als Todesart
. 
. *** (a) nur Variablen auf Individualebene
/************/
/* COLUMN 1 */ 
/************/
. regress suicide rel_protestant rel_otherchrist rel_jewish rel_islam rel_other rel_none rel_unknown german age age2 
> male famstand_unknown famstand_single famstand_widowed famstand_divorced 

      Source |       SS       df       MS              Number of obs =  885374
-------------+------------------------------           F( 14,885359) = 2104.43
       Model |  426.860318    14  30.4900227           Prob > F      =  0.0000
    Residual |  12827.5429885359  .014488521           R-squared     =  0.0322
-------------+------------------------------           Adj R-squared =  0.0322
       Total |  13254.4032885373  .014970417           Root MSE      =  .12037

------------------------------------------------------------------------------
     suicide |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
rel_protes~t |   .0017981     .00029     6.20   0.000     .0012296    .0023665
rel_otherc~t |  -.0002059   .0014819    -0.14   0.889    -.0031103    .0026985
  rel_jewish |   .0033212   .0055436     0.60   0.549    -.0075441    .0141865
   rel_islam |  (dropped)
   rel_other |  -.0104204   .0012669    -8.23   0.000    -.0129034   -.0079374
    rel_none |   .0044281   .0003923    11.29   0.000     .0036592     .005197
 rel_unknown |   .0031454   .0009567     3.29   0.001     .0012702    .0050205
      german |   .0042302   .0011976     3.53   0.000      .001883    .0065774
         age |  -.0031404   .0000475   -66.06   0.000    -.0032335   -.0030472
        age2 |   .0013997    .000037    37.82   0.000     .0013272    .0014723
        male |   .0071967   .0002898    24.83   0.000     .0066287    .0077647
famstand_u~n |  -.1333416   .0020092   -66.36   0.000    -.1372796   -.1294036
famstand_s~e |   .0088623   .0004947    17.91   0.000     .0078927    .0098319
famstand_w~d |    .007707   .0003422    22.52   0.000     .0070364    .0083777
famstand_d~d |   .0065111    .000568    11.46   0.000     .0053979    .0076243
       _cons |   .1542197   .0019149    80.54   0.000     .1504666    .1579728
------------------------------------------------------------------------------

. 
. * Kontrolliere fuer fixe Landkreiseffekte (Variante 1 via -absorb-)
/************/
/* COLUMN 2 */ 
/************/
. areg suicide rel_protestant rel_otherchrist rel_jewish rel_islam rel_other rel_none rel_unknown german age age2 mal
> e famstand_unknown famstand_single famstand_widowed famstand_divorced, absorb(ags)
(dropping rel_islam because it does not vary within category)

Linear regression, absorbing indicators                Number of obs =  885374
                                                       F( 14,884814) = 2101.69
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.0334
                                                       Adj R-squared =  0.0328
                                                       Root MSE      =  .12033

------------------------------------------------------------------------------
     suicide |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
rel_protes~t |   .0025201     .00034     7.41   0.000     .0018538    .0031864
rel_otherc~t |   .0003097   .0014899     0.21   0.835    -.0026105      .00323
  rel_jewish |   .0044137   .0055547     0.79   0.427    -.0064732    .0153007
   rel_islam |  (dropped)
   rel_other |  -.0113966    .001506    -7.57   0.000    -.0143483   -.0084448
    rel_none |   .0059679   .0004944    12.07   0.000     .0049989    .0069369
 rel_unknown |   .0043995   .0011017     3.99   0.000     .0022403    .0065587
      german |   .0041701   .0012186     3.42   0.001     .0017816    .0065586
         age |  -.0031224   .0000476   -65.60   0.000    -.0032157   -.0030291
        age2 |   .0013811   .0000371    37.26   0.000     .0013085    .0014538
        male |   .0071512   .0002902    24.64   0.000     .0065823      .00772
famstand_u~n |  -.1329624   .0020112   -66.11   0.000    -.1369044   -.1290204
famstand_s~e |   .0086981   .0004957    17.55   0.000     .0077265    .0096698
famstand_w~d |   .0078714   .0003424    22.99   0.000     .0072003    .0085425
famstand_d~d |    .006594   .0005699    11.57   0.000     .0054771    .0077109
       _cons |    .153375   .0019412    79.01   0.000     .1495703    .1571798
-------------+----------------------------------------------------------------
         ags |    F(545, 884814) =      1.995   0.000         (546 categories)

.
. * nur f�r Westdeutschland (da nur dort VZ 1987 verf�gbar)
/************/
/* COLUMN 5 */ 
/************/
. areg suicide rel_protestant rel_otherchrist rel_jewish rel_islam rel_other rel_none rel_unknown german age age2 mal
> e famstand_unknown famstand_single famstand_widowed famstand_divorced log_pop_tot_res_t share*, absorb(ags) vce(clu
> ster ags)
(dropping rel_islam because it does not vary within category)
(dropping log_pop_tot_res_t because it does not vary within category)
(dropping share_pop_under3_res_t because it does not vary within category)
(dropping share_pop_3_5_res_t because it does not vary within category)
(dropping share_pop_5_6_res_t because it does not vary within category)
(dropping share_pop_6_10_res_t because it does not vary within category)
(dropping share_pop_10_15_res_t because it does not vary within category)
(dropping share_pop_15_18_res_t because it does not vary within category)
(dropping share_pop_18_20_res_t because it does not vary within category)
(dropping share_pop_20_25_res_t because it does not vary within category)
(dropping share_pop_25_30_res_t because it does not vary within category)
(dropping share_pop_30_35_res_t because it does not vary within category)
(dropping share_pop_35_40_res_t because it does not vary within category)
(dropping share_pop_40_45_res_t because it does not vary within category)
(dropping share_pop_45_50_res_t because it does not vary within category)
(dropping share_pop_50_55_res_t because it does not vary within category)
(dropping share_pop_55_60_res_t because it does not vary within category)
(dropping share_pop_60_63_res_t because it does not vary within category)
(dropping share_pop_63_65_res_t because it does not vary within category)
(dropping share_pop_65_70_res_t because it does not vary within category)
(dropping share_pop_70_75_res_t because it does not vary within category)
(dropping share_pop_75andmore_res_t because it does not vary within category)
(dropping share_pop_tot_res_single_t because it does not vary within category)
(dropping share_pop_tot_res_married_t because it does not vary within category)
(dropping share_pop_tot_res_wid_t because it does not vary within category)
(dropping share_pop_tot_res_div_t because it does not vary within category)
(dropping share_foreigners because it does not vary within category)
(dropping share_pop_res_wrkmninc_1st_t because it does not vary within category)
(dropping share_pop_res_wrkmninc_2nd_t because it does not vary within category)
(dropping share_pop_res_wrkmninc_3rd_t because it does not vary within category)
(dropping share_pop_res_wrkmninc_oth_t because it does not vary within category)
(dropping share_pop_res_incsub_t because it does not vary within category)
(dropping share_pop_res_incpar_t because it does not vary within category)
(dropping share_school_prim because it does not vary within category)
(dropping share_school_lowsec because it does not vary within category)
(dropping share_school_uppsec because it does not vary within category)
(dropping share_voc_fach because it does not vary within category)
(dropping share_voc_uni because it does not vary within category)

Linear regression, absorbing indicators                Number of obs =  666261
                                                       F( 14,   326) =  266.95
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.0375
                                                       Adj R-squared =  0.0370
                                                       Root MSE      =  .11789

                                  (Std. Err. adjusted for 327 clusters in ags)
------------------------------------------------------------------------------
             |               Robust
     suicide |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
rel_protes~t |   .0021496    .000372     5.78   0.000     .0014178    .0028815
rel_otherc~t |   .0002909   .0013839     0.21   0.834    -.0024316    .0030134
  rel_jewish |   .0093053   .0077431     1.20   0.230    -.0059274     .024538
   rel_islam |  (dropped)
   rel_other |  -.0269936   .0035307    -7.65   0.000    -.0339394   -.0200477
    rel_none |   .0109761   .0012506     8.78   0.000     .0085159    .0134364
 rel_unknown |   .0184945   .0045606     4.06   0.000     .0095226    .0274665
      german |   .0016151   .0020585     0.78   0.433    -.0024345    .0056647
         age |  -.0033814    .000141   -23.99   0.000    -.0036587   -.0031041
        age2 |   .0015503   .0000948    16.35   0.000     .0013638    .0017368
        male |   .0070215   .0003581    19.61   0.000      .006317    .0077261
famstand_u~n |   -.141359   .0065179   -21.69   0.000    -.1541814   -.1285366
famstand_s~e |   .0112285   .0009184    12.23   0.000     .0094218    .0130352
famstand_w~d |   .0080391   .0003715    21.64   0.000     .0073082      .00877
famstand_d~d |    .007009   .0008609     8.14   0.000     .0053153    .0087027
log_pop_to~t |  (dropped)
sha~r3_res_t |  (dropped)
sha~_5_res_t |  (dropped)
shar~6_res_t |  (dropped)
sha~10_res_t |  (dropped)
sha~15_res_t |  (dropped)
shar~8_res_t |  (dropped)
sha~20_res_t |  (dropped)
sha~25_res_t |  (dropped)
sha~30_res_t |  (dropped)
sha~35_res_t |  (dropped)
sha~40_res_t |  (dropped)
sha~45_res_t |  (dropped)
sha~50_res_t |  (dropped)
sha~55_res_t |  (dropped)
sha~60_res_t |  (dropped)
sha~63_res_t |  (dropped)
sha~65_res_t |  (dropped)
sha~70_res_t |  (dropped)
sha~75_res_t |  (dropped)
shar~e_res_t |  (dropped)
share_po~e_t |  (dropped)
share_p~ed_t |  (dropped)
share_p~id_t |  (dropped)
share_po~v_t |  (dropped)
share_fore~s |  (dropped)
share_po~t_t |  (dropped)
share_p~nd_t |  (dropped)
share_p~rd_t |  (dropped)
share_po~h_t |  (dropped)
share_po~b_t |  (dropped)
share_po~r_t |  (dropped)
share_scho~m |  (dropped)
share_s~wsec |  (dropped)
share_s~psec |  (dropped)
share_voc_~h |  (dropped)
share_voc_~i |  (dropped)
       _cons |   .1653283   .0054315    30.44   0.000      .154643    .1760136
-------------+----------------------------------------------------------------
         ags |   absorbed                                     (327 categories)
		 
. *** (b) plus Variablen aus VZ 1987
/************/
/* COLUMN 6 */
/************/
. regress suicide rel_protestant rel_otherchrist rel_jewish rel_islam rel_other rel_none rel_unknown german age age2 
> male famstand_unknown famstand_single famstand_widowed famstand_divorced log_pop_tot_res_t share*, cluster(ags)

Linear regression                                      Number of obs =  666261
                                                       F( 47,   326) =       .
                                                       Prob > F      =       .
                                                       R-squared     =  0.0369
                                                       Root MSE      =  .11791

                                  (Std. Err. adjusted for 327 clusters in ags)
------------------------------------------------------------------------------
             |               Robust
     suicide |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
rel_protes~t |   .0021713   .0003408     6.37   0.000     .0015009    .0028417
rel_otherc~t |   .0003383   .0013798     0.25   0.806    -.0023761    .0030528
  rel_jewish |     .00836   .0077523     1.08   0.282    -.0068908    .0236108
   rel_islam |  (dropped)
   rel_other |  -.0271039   .0035347    -7.67   0.000    -.0340576   -.0201502
    rel_none |   .0106489   .0012114     8.79   0.000     .0082657    .0130321
 rel_unknown |   .0184102   .0045243     4.07   0.000     .0095097    .0273107
      german |   .0015453   .0020588     0.75   0.453    -.0025049    .0055956
         age |  -.0033827   .0001409   -24.01   0.000    -.0036598   -.0031056
        age2 |    .001551   .0000947    16.37   0.000     .0013646    .0017374
        male |   .0070224   .0003582    19.61   0.000     .0063178     .007727
famstand_u~n |  -.1413783   .0065123   -21.71   0.000    -.1541897   -.1285668
famstand_s~e |   .0112121   .0009175    12.22   0.000     .0094071     .013017
famstand_w~d |   .0080274   .0003707    21.66   0.000     .0072982    .0087565
famstand_d~d |   .0070074   .0008601     8.15   0.000     .0053154    .0086994
log_pop_to~t |  -.0005275   .0009256    -0.57   0.569    -.0023485    .0012934
sha~r3_res_t |   1.106834   .6298338     1.76   0.080    -.1322174    2.345886
sha~_5_res_t |  (dropped)
shar~6_res_t |   2.687555    .943727     2.85   0.005     .8309916    4.544118
sha~10_res_t |   .9645112   .5357589     1.80   0.073    -.0894698    2.018492
sha~15_res_t |   .5492675   .4890085     1.12   0.262    -.4127429    1.511278
shar~8_res_t |    1.31777   .4916865     2.68   0.008     .3504911    2.285049
sha~20_res_t |   .6560224   .5024307     1.31   0.193    -.3323932    1.644438
sha~25_res_t |   .9572602   .4385949     2.18   0.030     .0944266    1.820094
sha~30_res_t |   .8865327   .4480884     1.98   0.049      .005023    1.768042
sha~35_res_t |   1.064604   .4728252     2.25   0.025     .1344299    1.994777
sha~40_res_t |   .6723428   .4405848     1.53   0.128    -.1944054    1.539091
sha~45_res_t |    .847552    .439253     1.93   0.055    -.0165761     1.71168
sha~50_res_t |   1.099842   .4664727     2.36   0.019     .1821654    2.017519
sha~55_res_t |   .8028413   .4447359     1.81   0.072    -.0720731    1.677756
sha~60_res_t |   .9805729   .4968209     1.97   0.049     .0031933    1.957952
sha~63_res_t |   .8712311   .4864739     1.79   0.074    -.0857932    1.828255
sha~65_res_t |   .5509281   .5258931     1.05   0.296    -.4836443    1.585501
sha~70_res_t |    1.29907   .4986202     2.61   0.010     .3181504    2.279989
sha~75_res_t |   .9503095    .439234     2.16   0.031     .0862188      1.8144
shar~e_res_t |   .9474508    .459571     2.06   0.040     .0433517     1.85155
share_po~e_t |  -.1023426   .0731282    -1.40   0.163    -.2462053    .0415201
share_p~ed_t |  -.1445334   .0635306    -2.28   0.024    -.2695151   -.0195518
share_p~id_t |  -.1733315    .113894    -1.52   0.129    -.3973915    .0507284
share_po~v_t |  (dropped)
share_fore~s |  -.0270688   .0143037    -1.89   0.059     -.055208    .0010704
share_po~t_t |  (dropped)
share_p~nd_t |  -.0353077   .0327576    -1.08   0.282    -.0997508    .0291354
share_p~rd_t |  -.0293931   .0399445    -0.74   0.462    -.1079747    .0491884
share_po~h_t |  -.0399528   .0372978    -1.07   0.285    -.1133276    .0334219
share_po~b_t |  -.1227869   .0395969    -3.10   0.002    -.2006847   -.0448891
share_po~r_t |  -.0612559   .0350011    -1.75   0.081    -.1301124    .0076005
share_scho~m |   6.51e-11   1.51e-08     0.00   0.997    -2.96e-08    2.97e-08
share_s~wsec |  -3.68e-08   5.82e-08    -0.63   0.528    -1.51e-07    7.78e-08
share_s~psec |  -6.42e-08   9.96e-08    -0.64   0.520    -2.60e-07    1.32e-07
share_voc_~h |   7.41e-08   9.44e-08     0.79   0.433    -1.12e-07    2.60e-07
share_voc_~i |   1.15e-07   1.53e-07     0.75   0.454    -1.86e-07    4.15e-07
       _cons |  -.5609905   .4460257    -1.26   0.209    -1.438442    .3164614
------------------------------------------------------------------------------

. * (nur) Protestantenanteil statt individueller Religionszugehoerigkeit
/************/
/* COLUMN 7 */ 
/************/
. regress suicide prot_in_protcath rel_otherchrist rel_jewish rel_islam rel_other rel_none rel_unknown german age age
> 2 male famstand_unknown famstand_single famstand_widowed famstand_divorced log_pop_tot_res_t share*, cluster(ags)

Linear regression                                      Number of obs =  666261
                                                       F( 47,   326) =       .
                                                       Prob > F      =       .
                                                       R-squared     =  0.0368
                                                       Root MSE      =  .11791

                                  (Std. Err. adjusted for 327 clusters in ags)
------------------------------------------------------------------------------
             |               Robust
     suicide |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
prot_in_pr~h |   .0024375   .0011992     2.03   0.043     .0000784    .0047967
rel_otherc~t |  -.0009283   .0013116    -0.71   0.480    -.0035087     .001652
  rel_jewish |   .0074173    .007707     0.96   0.337    -.0077445    .0225791
   rel_islam |  (dropped)
   rel_other |   -.027874   .0035209    -7.92   0.000    -.0348007   -.0209474
    rel_none |    .009269   .0012858     7.21   0.000     .0067395    .0117986
 rel_unknown |   .0172062    .004574     3.76   0.000     .0082078    .0262046
      german |   .0019898   .0020488     0.97   0.332    -.0020407    .0060203
         age |  -.0033874    .000141   -24.03   0.000    -.0036648   -.0031101
        age2 |   .0015556   .0000948    16.41   0.000      .001369    .0017421
        male |   .0070281   .0003586    19.60   0.000     .0063227    .0077335
famstand_u~n |  -.1414344   .0065153   -21.71   0.000    -.1542517   -.1286171
famstand_s~e |   .0111377   .0009167    12.15   0.000     .0093343     .012941
famstand_w~d |   .0080313   .0003707    21.67   0.000      .007302    .0087606
famstand_d~d |    .007064   .0008612     8.20   0.000     .0053698    .0087583
log_pop_to~t |  -.0005098   .0009278    -0.55   0.583    -.0023351    .0013155
sha~r3_res_t |   1.113175   .6304221     1.77   0.078    -.1270338    2.353384
sha~_5_res_t |  (dropped)
shar~6_res_t |   2.722762   .9601353     2.84   0.005     .8339194    4.611606
sha~10_res_t |   .9641393   .5365186     1.80   0.073    -.0913364    2.019615
sha~15_res_t |   .5458042   .4881673     1.12   0.264    -.4145514     1.50616
shar~8_res_t |   1.299613   .4911244     2.65   0.009     .3334399    2.265786
sha~20_res_t |   .6631235   .5028595     1.32   0.188    -.3261356    1.652383
sha~25_res_t |    .954714   .4387568     2.18   0.030      .091562    1.817866
sha~30_res_t |   .8804133   .4477745     1.97   0.050    -.0004789    1.761305
sha~35_res_t |   1.064498   .4728437     2.25   0.025     .1342878    1.994708
sha~40_res_t |   .6649234   .4401674     1.51   0.132    -.2010036     1.53085
sha~45_res_t |   .8515632   .4409112     1.93   0.054    -.0158271    1.718953
sha~50_res_t |   1.091953   .4661314     2.34   0.020     .1749479    2.008958
sha~55_res_t |   .8060574   .4462123     1.81   0.072    -.0717617    1.683876
sha~60_res_t |   .9825298   .4970131     1.98   0.049     .0047721    1.960287
sha~63_res_t |   .8825662   .4941095     1.79   0.075    -.0894793    1.854612
sha~65_res_t |   .5769632   .5305441     1.09   0.278    -.4667589    1.620685
sha~70_res_t |   1.283212   .4979517     2.58   0.010     .3036078    2.262816
sha~75_res_t |   .9723589   .4451469     2.18   0.030      .096636    1.848082
shar~e_res_t |   .9388247   .4593991     2.04   0.042     .0350638    1.842586
share_po~e_t |  -.1009348   .0734835    -1.37   0.171    -.2454965    .0436269
share_p~ed_t |  -.1471835   .0638966    -2.30   0.022    -.2728852   -.0214818
share_p~id_t |  -.1714667   .1139276    -1.51   0.133    -.3955928    .0526593
share_po~v_t |  (dropped)
share_fore~s |  -.0266802   .0143208    -1.86   0.063    -.0548531    .0014927
share_po~t_t |  (dropped)
share_p~nd_t |   -.033187   .0324521    -1.02   0.307    -.0970289     .030655
share_p~rd_t |  -.0262421   .0397617    -0.66   0.510    -.1044641    .0519798
share_po~h_t |  -.0377865   .0370698    -1.02   0.309    -.1107128    .0351397
share_po~b_t |  -.1245275   .0402181    -3.10   0.002    -.2036473   -.0454078
share_po~r_t |  -.0593173     .03474    -1.71   0.089    -.1276601    .0090256
share_scho~m |   9.20e-10   1.54e-08     0.06   0.952    -2.94e-08    3.12e-08
share_s~wsec |  -4.04e-08   5.96e-08    -0.68   0.498    -1.58e-07    7.67e-08
share_s~psec |  -6.18e-08   1.00e-07    -0.62   0.538    -2.59e-07    1.36e-07
share_voc_~h |   7.28e-08   9.48e-08     0.77   0.443    -1.14e-07    2.59e-07
share_voc_~i |   1.14e-07   1.53e-07     0.74   0.458    -1.87e-07    4.15e-07
       _cons |  -.5612644   .4461689    -1.26   0.209    -1.438998    .3164692
------------------------------------------------------------------------------


/************/
/* COLUMN 8 */ 
/************/
. regress suicide rel_protestant prot_in_protcath rel_otherchrist rel_jewish rel_islam rel_other rel_none rel_unknown
>  german age age2 male famstand_unknown famstand_single famstand_widowed famstand_divorced log_pop_tot_res_t share*,
>  cluster(ags)

Linear regression                                      Number of obs =  666261
                                                       F( 48,   326) =       .
                                                       Prob > F      =       .
                                                       R-squared     =  0.0369
                                                       Root MSE      =  .11791

                                  (Std. Err. adjusted for 327 clusters in ags)
------------------------------------------------------------------------------
             |               Robust
     suicide |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
rel_protes~t |   .0021318   .0003694     5.77   0.000     .0014051    .0028586
prot_in_pr~h |   .0004617   .0012953     0.36   0.722    -.0020865      .00301
rel_otherc~t |   .0003051   .0013837     0.22   0.826     -.002417    .0030272
  rel_jewish |   .0083497   .0077542     1.08   0.282     -.006905    .0236043
   rel_islam |  (dropped)
   rel_other |   -.027122   .0035335    -7.68   0.000    -.0340733   -.0201707
    rel_none |   .0106198   .0012173     8.72   0.000     .0082251    .0130146
 rel_unknown |   .0183843   .0045305     4.06   0.000     .0094717     .027297
      german |   .0015501   .0020604     0.75   0.452    -.0025033    .0056035
         age |  -.0033826   .0001409   -24.01   0.000    -.0036598   -.0031055
        age2 |    .001551   .0000947    16.37   0.000     .0013646    .0017374
        male |   .0070224   .0003582    19.61   0.000     .0063178     .007727
famstand_u~n |  -.1413731   .0065126   -21.71   0.000    -.1541852   -.1285611
famstand_s~e |   .0112135   .0009172    12.23   0.000     .0094092    .0130178
famstand_w~d |   .0080278   .0003706    21.66   0.000     .0072987    .0087568
famstand_d~d |   .0070092   .0008603     8.15   0.000     .0053166    .0087017
log_pop_to~t |   -.000528   .0009268    -0.57   0.569    -.0023513    .0012953
sha~r3_res_t |   1.115311   .6310438     1.77   0.078    -.1261211    2.356743
sha~_5_res_t |  (dropped)
shar~6_res_t |   2.734391   .9600749     2.85   0.005     .8456669    4.623115
sha~10_res_t |   .9665083   .5369611     1.80   0.073    -.0898379    2.022855
sha~15_res_t |   .5505906   .4889699     1.13   0.261    -.4113441    1.512525
shar~8_res_t |   1.302186   .4910148     2.65   0.008      .336229    2.268144
sha~20_res_t |   .6638325   .5036786     1.32   0.188    -.3270381    1.654703
sha~25_res_t |   .9597475   .4389486     2.19   0.029     .0962181    1.823277
sha~30_res_t |   .8838452    .448076     1.97   0.049     .0023599     1.76533
sha~35_res_t |    1.06906    .473295     2.26   0.025     .1379619    2.000158
sha~40_res_t |   .6688187   .4406774     1.52   0.130    -.1981115    1.535749
sha~45_res_t |   .8578002   .4412031     1.94   0.053    -.0101644    1.725765
sha~50_res_t |   1.095661   .4664922     2.35   0.019     .1779459    2.013376
sha~55_res_t |   .8099736   .4466896     1.81   0.071    -.0687844    1.688732
sha~60_res_t |   .9831752   .4971157     1.98   0.049     .0052156    1.961135
sha~63_res_t |   .8873704   .4941549     1.80   0.073    -.0847644    1.859505
sha~65_res_t |   .5775455   .5307146     1.09   0.277    -.4665121    1.621603
sha~70_res_t |   1.290045    .498438     2.59   0.010     .3094844    2.270606
sha~75_res_t |   .9763975   .4456902     2.19   0.029     .0996055    1.853189
shar~e_res_t |   .9399529    .459666     2.04   0.042     .0356668    1.844239
share_po~e_t |  -.1020131   .0733408    -1.39   0.165    -.2462941    .0422679
share_p~ed_t |  -.1477999   .0637733    -2.32   0.021     -.273259   -.0223408
share_p~id_t |  -.1723185   .1140865    -1.51   0.132    -.3967572    .0521202
share_po~v_t |  (dropped)
share_fore~s |  -.0279125   .0143606    -1.94   0.053    -.0561636    .0003387
share_po~t_t |  (dropped)
share_p~nd_t |  -.0341515   .0324491    -1.05   0.293    -.0979876    .0296845
share_p~rd_t |  -.0271948   .0397227    -0.68   0.494      -.10534    .0509504
share_po~h_t |  -.0393249   .0370919    -1.06   0.290    -.1122945    .0336447
share_po~b_t |    -.12485   .0401896    -3.11   0.002    -.2039137   -.0457862
share_po~r_t |  -.0600869   .0347185    -1.73   0.084    -.1283875    .0082136
share_scho~m |   6.97e-10   1.54e-08     0.05   0.964    -2.95e-08    3.09e-08
share_s~wsec |  -4.02e-08   5.95e-08    -0.68   0.500    -1.57e-07    7.69e-08
share_s~psec |  -5.93e-08   1.00e-07    -0.59   0.555    -2.56e-07    1.38e-07
share_voc_~h |   7.44e-08   9.49e-08     0.78   0.434    -1.12e-07    2.61e-07
share_voc_~i |   1.10e-07   1.53e-07     0.72   0.475    -1.92e-07    4.11e-07
       _cons |  -.5627912   .4466577    -1.26   0.209    -1.441486     .315904
------------------------------------------------------------------------------


. /*------------------------------------------------------------------------------------------------------*/
. /*------------------------------------------------------------------------------------------------------*/
. /*------------------------------------------------------------------------------------------------------*/
. 
. /********************/
. /* 1998-2009: ICD10 */
. /********************/
. 
. **** Datensatz einlesen
. use "$root/KDFV/tus080_2009.dta", clear

. 
. *** Generiere Amtlichen Gemeindeschluessel Wohnort der/des Verstorbenen numerisch statt string
. gen ags=real(ef304)

. 
. *** Generiere Dummyvariable fuer Maenner
. gen male=(ef306==1)

. 
. *** Generiere Dummyvariablen fuer Familienstand
. gen famstand_unknown  = (ef310==0 | ef310==9) 

. gen famstand_single   = (ef310==1)

. gen famstand_married  = (ef310==2) /* reference category */

. gen famstand_widowed  = (ef310==3)

. gen famstand_divorced = (ef310==4)

. 
. *** Generiere Dummyvariablen fuer Religionszugehoerigkeit
. gen rel_unknown = (ef313==0) 

. gen rel_protestant = (ef313==1 | ef313==2)

. gen rel_catholic = (ef313==3 | ef313==4 | ef313==5) /* reference category */

. gen rel_otherchrist = (ef313==6)

. gen rel_jewish = (ef313==7)

. gen rel_islam = (ef313==8)

. gen rel_other = (ef313==9 | ef313==10)

. gen rel_none  = (ef313==11)

. gen rel_pro_no = rel_protestant + rel_none

. 
. *** Protestantenanteil an "Protestanten plus Katholiken" auf Landkreisebene aus TUS
. bysort ags: egen n_protestant_ags = sum(rel_protestant)

. bysort ags: egen n_catholic_ags = sum(rel_catholic)

. bysort ags: gen prot_in_protcath_ags = n_protestant_ags/(n_catholic_ags+n_protestant_ags)
(* missing values generated)

. 
. *** Deutsche vs Auslaender
. gen german = (ef314==0)

. 
. **** Generiere Selbstmord-Dummy
. *** Selbstmord-Codes in ICD10 sind Codes X60, X61, ..., X84, d.h. 25 Codes */
. gen suicide=0

. foreach num of numlist 60/84 {
  2.   replace suicide=1 if ef317=="X`num'"
  3. }
(* real changes made)
(* real changes made)
(* real changes made)
(* real changes made)
(* real changes made)
(* real changes made)
(* real change made)
(* real changes made)
(* real changes made)
(* real changes made)
(* real changes made)
(* real changes made)
(* real changes made)
(* real changes made)
(* real changes made)
(* real changes made)
(* real changes made)
(* real change made)
(* real changes made)
(* real changes made)
(* real changes made)
(* real changes made)
(* real changes made)
(* real changes made)
(* real changes made)

. 
. *** Generiere alter^2 auf Basis von Alter: Sterbealter - in Jahren
. rename ef320u3 age

. gen age2 = age * age / 100

. 
. *** Sortiere nach AGS
. sort ags

. 
. *** Verknuepfe mit Kreisdaten der VZ 1987
. merge ags using "$root/KDFV/VZ1987.dta"
variable ags does not uniquely identify observations in the master data

. 
. *** Benenne Variable Protestantenanteil um
. ren share_prot_in_protcath prot_in_protcath

. 
. *** Betrachte Korrelation zwischen Protestantenanteil in TUS und Protestantenanteil in VZ1987
. pwcorr prot_in_protcath prot_in_protcath_ags, sig

             | prot_i~h prot_i~s
-------------+------------------
prot_in_pr~h |   1.0000 
             |
             |
prot_in_pr~s |   0.9973   1.0000 
             |   0.0000
             |

. 
. *** Generiere log(Bevoelkerungszahl) im Landkreis
. gen log_pop_tot_res_t = log(pop_tot_res_t)
(* missing values generated)

. 
. *** Tabelliere Religionszugehoerigkeit der Verstorben insgesamt, sowie fuer die Selbstmoerder
. *tab ef313
. *tab ef313 if suicide==1
. 
. 
. 
. *** Regressionen: was "erklaert" Selbstmord als Todesart
. 
. *** (a) nur Variablen auf Individualebene
/************/
/* COLUMN 3 */ 
/************/
. regress suicide rel_protestant rel_otherchrist rel_jewish rel_islam rel_other rel_none rel_unknown german age age2 
> male famstand_unknown famstand_single famstand_widowed famstand_divorced 

      Source |       SS       df       MS              Number of obs =  854544
-------------+------------------------------           F( 14,854529) = 1793.46
       Model |  271.559066    14  19.3970761           Prob > F      =  0.0000
    Residual |  9242.09908854529  .010815431           R-squared     =  0.0285
-------------+------------------------------           Adj R-squared =  0.0285
       Total |  9513.65814854543  .011133036           Root MSE      =    .104

------------------------------------------------------------------------------
     suicide |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
rel_protes~t |   .0002655   .0002644     1.00   0.315    -.0002527    .0007837
rel_otherc~t |  -.0006577    .001426    -0.46   0.645    -.0034525    .0021371
  rel_jewish |   .0040461   .0029513     1.37   0.170    -.0017383    .0098305
   rel_islam |  -.0108808   .0017183    -6.33   0.000    -.0142487    -.007513
   rel_other |   .0061249   .0063814     0.96   0.337    -.0063824    .0186323
    rel_none |   .0016101   .0003158     5.10   0.000     .0009912     .002229
 rel_unknown |  -.0021039   .0006089    -3.46   0.001    -.0032974   -.0009104
      german |   .0011208   .0008018     1.40   0.162    -.0004507    .0026924
         age |  -.0026508   .0000421   -62.95   0.000    -.0027334   -.0025683
        age2 |   .0011326   .0000313    36.18   0.000     .0010712    .0011939
        male |   .0073762   .0002522    29.24   0.000     .0068818    .0078706
famstand_u~n |  (dropped)
famstand_s~e |   .0074022   .0004255    17.40   0.000     .0065683    .0082361
famstand_w~d |   .0060977   .0002979    20.47   0.000     .0055138    .0066815
famstand_d~d |   .0042906   .0004336     9.90   0.000     .0034409    .0051404
       _cons |   .1371392    .001664    82.41   0.000     .1338777    .1404006
------------------------------------------------------------------------------

. 
. * Kontrolliere fuer fixe Landkreiseffekte (Variante 1 via -absorb-)
/************/
/* COLUMN 4 */ 
/************/
. areg suicide rel_protestant rel_otherchrist rel_jewish rel_islam rel_other rel_none rel_unknown german age age2 mal
> e famstand_unknown famstand_single famstand_widowed famstand_divorced, absorb(ags)
(dropping famstand_unknown because it does not vary within category)

Linear regression, absorbing indicators                Number of obs =  854544
                                                       F( 14,854007) = 1805.64
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.0300
                                                       Adj R-squared =  0.0294
                                                       Root MSE      =  .10395

------------------------------------------------------------------------------
     suicide |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
rel_protes~t |   .0014313   .0003013     4.75   0.000     .0008407    .0020219
rel_otherc~t |   -.000207   .0014316    -0.14   0.885    -.0030128    .0025988
  rel_jewish |   .0052171   .0029588     1.76   0.078    -.0005821    .0110162
   rel_islam |   -.009749    .001723    -5.66   0.000    -.0131261   -.0063719
   rel_other |   .0072146   .0063822     1.13   0.258    -.0052943    .0197235
    rel_none |   .0038079   .0003957     9.62   0.000     .0030325    .0045834
 rel_unknown |   .0002764   .0006539     0.42   0.673    -.0010052    .0015581
      german |   .0016339   .0008061     2.03   0.043     .0000539    .0032139
         age |  -.0026274   .0000422   -62.33   0.000      -.00271   -.0025448
        age2 |   .0011122   .0000313    35.49   0.000     .0010508    .0011736
        male |   .0073146   .0002525    28.97   0.000     .0068198    .0078094
famstand_u~n |  (dropped)
famstand_s~e |   .0072503   .0004265    17.00   0.000     .0064145    .0080862
famstand_w~d |   .0062357    .000298    20.92   0.000     .0056516    .0068199
famstand_d~d |   .0043579   .0004349    10.02   0.000     .0035056    .0052103
       _cons |   .1350266   .0016744    80.64   0.000     .1317448    .1383084
-------------+----------------------------------------------------------------
         ags |    F(522, 854007) =      2.446   0.000         (523 categories)

. 
. log close
       log:  D:/data/23211 Todesursachenstatistik/Nutzungen/1813-2011 ifo Muenchen - Knoche//Analysen/20121004/201210
> 04_tus-analyse.log
  log type:  text
 closed on:   4 Oct 2012, 22:22:19
---------------------------------------------------------------------------------------------------------------------
